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O Abstract 

(N ■ 

Digital radio antenna arrays, like LOPES (L0FAR PrototypE Station), detect high-energy cosmic rays via the radio emission 
from atmospheric extensive air showers. LOPES is an array of dipole antennas placed within and triggered by the KASCADE- 

ry | Grande experiment on site of the Karlsruhe Institute of Technology, Germany. The antennas are digitally combined to build a 
radio interferometer by forming a beam into the air shower arrival direction which allows measurements even at low signal-to-noise 
ratios in individual antennas. This technique requires a precise time calibration. A combination of several calibration steps is used 
to achieve the necessary timing accuracy of about 1 ns. The group delays of the setup are measured, the frequency dependence 

|* ""of these delays (dispersion) is corrected in the subsequent data analysis, and variations of the delays with time are monitored. 

J> We use a transmitting reference antenna, a beacon, which continuously emits sine waves at known frequencies. Variations of the 

HH relative delays between the antennas can be detected and corrected for at each recorded event by measuring the phases at the beacon 

r £^ frequencies. 
Oh' 
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1. Introduction 

For the study of ultra-high energy particles from the cosmos 
the measurement of the radio emission from secondary particle 
showers generated in air or dense media is evolving as a new 
technique [1]. First measurements of the radio emission of cos- 
mic ray air showers had been done already in the 1960's [2], but 
with the analog electronics available at that time, the technique 
could not be competitive with traditional methods like the de- 
tection of secondary particles on ground or the measurement of 
fluorescence light emitted by air showers. Recently, the radio 
detection method experienced a revival because of the avail- 
ability of fast digital electronics. Pioneering experiments like 
LOPES H and CODALEMA 0] have proven that radio detec- 
tion of cosmic ray air showers is possible with modern, digital 
antenna arrays. Due to the short duration of typically less than 
100 ns of the air shower induced radio pulse, the experimen- 
tal procedures are significantly different from those of classical 
radio astronomy. 

The main goal of the investigations is the detailed under- 
standing of the shower radio emission and the correlation of 
the measured field strengths with the primary cosmic ray char- 
acteristics. The sensitivity of the measurements to the direction 
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of the shower axis, the energy and mass of the primary particle 
are of particular interest. Radio antenna arrays can derive the 
energy of the primary particle by measuring the amplitude of 
the field strength, and reconstruct the direction of the incoming 
primary particle by measuring pulse arrival times - with the re- 
markable difference to other distributed sensor networks, that 
with LOPES, the arrival direction is reconstructed using digital 
interferometry which demands a precise time calibration. An- 
other goal is the optimization of the hardware (antenna design 
and electronics) for a large scale application of the detection 
technique including a self-trigger mechanism for stand-alone 
radio operation Q 01 • 

LOPES was built as a prototype station of the astronomical 
radio telescope LOFAR |7l |8[] aiming to investigate the new 
detection method in detail. LOPES is a phased array of ra- 
dio antennas. Featuring a precise time calibration, it can be 
used for interferometric measurements, e.g. when forming a 
cross-correlation beam into the air shower direction [9]. Thus, 
LOPES is sensitive to the coherence of the radio signal emitted 
by air showers, allowing to perform measurements even at low 
signal-to-noise ratios in individual antennas. 

This paper describes methods for the calibration and con- 
tinuous monitoring of the timing of a radio antenna array like 
LOPES and shows that it is possible to achieve a timing ac- 
curacy in the order of 1 ns by combining these methods for 
such kind of arrays. Beside the measurement and correction of 
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group delays and frequency dependent dispersion of the setup, 
we use a transmitting reference antenna, a beacon, which con- 
tinuously emits sine waves at known frequencies. This way, 
variations of the relative delays between the antennas can be 
detected and corrected for in the subsequent analysis of each 
recorded event by measuring the relative phases at the beacon 
frequencies. This is different from the time calibration in other 
experiments, like ANTARES d, ANITA HI and AURA HI 
which determine the arrival times of pulses emitted by a beacon. 
In addition, AURA has the capability to measure frequency 
shifts of constant waves for calibration |13]. The use of phase 
differences of a continuously emitting beacons is reported for 
ionospheric TEC measurements 11411 . where the measurement 
of phases of a beacon signal is used for atmospheric monitor- 
ing, not for time calibration. Where the individual methods de- 
scribed in this work are more or less standard in sensor based 
experiments, their combination to achieve the possibility of in- 
terferometric measurements is new and applied for the first time 
in LOPES. 

2. The LOPES antenna array 

The main component of LOPES consists of 30 amplitude cal- 
ibrated, inverted V-shape dipole antennas lfl5l [l6ll . The anten- 
nas are placed in co-location with the particle air shower ex- 
periment KASCADE-Grande d (fig. ffl). KASCADE- 
Grande consists mainly of stations equipped with scintilla- 
tion detectors on an area of 700 x 700 m 2 , where 252 sta- 
tions compose the KASCADE array, and further 37 large sta- 
tions the Grande array. Besides the 30 LOFAR-type antennas, 
LOPES consists also of newly designed antennas forming the 
LOPES STAR array [19]. The main purpose of LOPES STAR is 
to optimize the hardware for an application of this measuring 
technique to large scales, e.g. at the Pierre Auger Observa- 
tory 12011 . All antennas are optimized to measure in the range 
of 40 to 80 MHz which is less polluted by strong interference 
than, e.g. the FM band. The positions of the antennas have been 
determined by differential GPS measurements with a relative 
accuracy of a few cm. 

Whenever KASCADE-Grande measures a high-energy 
event, a trigger signal is send to LOPES which then stores the 
digitally recorded radio signal as a trace of 2 16 samples with 
a sampling frequency of 80 MHz, where the trigger time is 
roughly in the middle of the trace. As a band-pass filter is 
used to restrict the frequency band to 40 to 80 MHz, LOPES 
is operating in the second Nyquist domain and contains the 
complete information of the radio signal within this frequency 
band. Recovery of the full information is possible by an up- 
sampling procedure, i.e. the correct interpolation between the 
sam pled data points which is done by a zero-padding algorithm 
1 21), |22j, |23J . This way, sample spacings of 0.1 ns can be ob- 
tained within reasonable computing time, which is consider- 
ably smaller than the uncertainties of the timing introduced by 
other sources (see below). Thus, the sampling rate does not 
contribute significantly to systematic uncertainties. 

More details of the experimental set-up, the amplitude cali- 
bration, the operation, and the analysis procedures of LOPES 
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Figure 1 : Current setup of the LOPES experiment which is co-located with the 
KASCADE-Grande experiment at Karlsruhe Institute of Technology, Germany. 
Upward triangles mark east-west oriented antennas, downward triangles north- 
south oriented antennas, respectively. A star indicates an east-west oriented and 
a north-south oriented antenna at the same place. 



can be found in references, e.g. fl fllill . 



3. Need for a precise time calibration 

The angular resolution, respectively source location, of LOPES 
is limited to about 1° ll24Tl due to the uncertainties of the emis- 
sion mechanism of the radio pulse, and thus, by the uncertain- 
ties in the shape of the wave front of the radio emission. Con- 
sequently, for LOPES, improving the accuracy of the time cal- 
ibration to about 1 ns is not expected to significantly improve 
the angular resolution. Instead, this good timing resolution is a 
necessary requirement to enable the use of LOPES as a digital 
radio interferometer. Hence, this is the most important among 
several reasons why a precise time calibration with a relative 
accuracy in the order of or below 1 ns is desirable for a radio air 
shower array: 

• Interferometry: A timing precision which is at least an or- 
der of magnitude better than the period of the filter ringing 
(~ 17 ns for LOPES) allows one to perform interferomet- 
ric measurements if the baselines of the interferometer are 
adequate for the angular scale of the observed source. As 
the distance of the source of radio emission from cosmic 
ray air showers to the LOPES antenna array (several km) 
is much larger than the extension of the source region and 
the lateral extension of the array (~200m), the angular 
extension of the source is small. Hence, one expects that 
every antenna detects the same radio pulse just at a dif- 
ferent time. Thus, LOPES should see coherent radio sig- 
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Figure 2: Calibrated, up-sampled traces of an example event, after correction 
for geometrical delays. The radio pulse induced by a cosmic ray air shower 
can clearly be distinguished from the noise, as it is (in contrast to the noise) 
coherently detected in all antennas. 
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Figure 3: Influence of an additional timing uncertainty on the height of the 
cross-correlation beam of the example event (figure [2): The x-axis shows the 
width of the Gaussian distribution of the additional timing uncertainty added to 
each antenna. The error bars are the RMS of 100 repetitions which have been 
performed for each uncertainty. 



nals from air showers on the ground, which has been ex- 
perimentally verified H, and can be expemplarily seen in 
figure[2] This coherence is measurable, e.g., by forming a 
cross-correlation beam into the air shower direction [25], 
and can be used to distinguish between noise (e.g. thermal 
noise and noise originating from the KASCADE particle 
detectors) and air shower signals. 

The requirement of a timing precision in the order of 1 ns 
for the interferometric cross-correlation beam analysis, 
can be quantitatively verified by adding an additional and 
random timing uncertainty to each antenna, and studying 
the influence on the reconstructed cross-correlation beam 
which is a measure for the coherence. This has been done 
for the example event (figure |2]i by shifting the traces of 
each antenna by an additional time taken from a Gaussian 
random distribution (see figure[3]l. The height of the cross- 
correlation beam decreases significantly when the added 
uncertainty is larger than 1 ns. For uncertainties > 5 ns 
the height is not reduced further, as the analysis always 
finds a random correlation between some antennas. As 
most of the LOPES events are closer to the noise than the 
shown example, reconstructing the cross-correlation beam 
correctly is important, because a reduced height can lead 
to a signal-to-noise ratio below the detection threshold. 

Polarization studies: Different models for the radio emis- 
sion of air showers can, among others, also be tested by 
their predictions on the polarization of the radio signal 
(e.g., the geo-synchrotron model |26l 27 ] predicts predom- 



inantly linear polarization of the electric field in a direction 
depending on the geometry of the air shower [28]). The ca- 
pability of any antenna array to reconstruct the time depen- 
dence of the polarization vector at each antenna position, 
and thus, to distinguish between linearly and circularly po- 
larized signals, depends strongly on the relative timing ac- 
curacy between the different polarization channels of each 



antenna. 



• Lateral distribution of arrival times: According to simu- 
lations, the lateral distribution of the pulse arrival times 
should contain information about the mass of the primary 
cosmic ray particle 12911 . Only a precise relative timing, 
even between distant antennas (~200m for LOPES), can 
enable us to reveal this information, and to measure the 
shape of the radio wave front in detail. 



As stable clocks for the DAQ electronics and the trigger signal 
of LOPES are distributed via cables, the time calibration is ba- 
sically reduced to the measurement of the electronics and cable 
group delays, their dependence on the frequency (dispersion), 
and their variations with time. Originally, the delays were mea- 
sured with the radio emission from solar burst events, and their 
variations were monitored by measuring the phase of the carri- 
ers of a television transmitter yl Pol . 

Meanwhile, we have developed new methods for the time 
calibration which do not depend on external sources out of our 
control. Namely, we measure the delays with a reference pulse 
emitted at a known time, correct for the dispersion of the ana- 
log electronics and have set up an emitting antenna (beacon) 
which continuously transmits two narrow band reference sig- 
nals to monitor variations of the delays with time. 

These three methods for calibration and monitoring of the 
timing are combined to achieve a timing accuracy in the or- 
der of 1 ns for each event measured with LOPES. Nevertheless, 
these methods are in principle independent from each other, and 
for other experiments one might, e.g., determine delays by an- 
other method, but still use the beacon method to continuously 
monitor the relative timing. 
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Figure 4: Setup for the delay measurements: the coaxial cable from the LOPES 
antenna is disconnected from the antenna and connected to a cable from the 
pulse generator instead which is triggered simultaneously with the DAQ elec- 
tronics. 



4. Delay measurements 

For LOPES, as a digital radio interferometer, mainly the relative 
timing between the different antennas is of importance, and the 
absolute event time has to be known only roughly to combine 
the LOPES events with the corresponding KASCADE-Grande 
events. Thus, the determination of the pulse arrival times at 
each antenna, and therefore the measurement of the delays, is 
most important on a relative basis. Hereby, the delay of each 
channel (antenna and its analog electronics) is different, e.g., 
because different cable lengths are used. 

We define the absolute delay t of a channel as the time be- 
tween the arrival time frj of a radio pulse at an antenna and 
the time f t when it appears in the digitally measured trace: 
r = f t - to. The more important relative delay AT m n between 
two antennas m and n is the difference between the absolute 
delays of these antennas: AT m n = r m - r n . 

Using solar bursts all relative delays Ar m n could be deter- 
mined directly. Measuring the delays with respect to a common 
reference time f re f is equivalent if the difference t re f -to = const 
is the same for all antennas. These delays f measured with re- 
spect to f le f are related to the absolute delays by f = t t - f le f = 
t - (h-ef - to), and the relative delays can be easily derived from 
the measured delays f by Ar m n = r m - r n = f m — f n . 

For each antenna the measurement of the delay f is per- 
formed as follows: We disconnect the cable from the antenna 
and connect it to a pulse generator instead, which emits a short 
calibration pulse at a fixed time after a normal KASCADE- 
Grande trigger (fig. @}. As reference time f re f we define the 
zero point of the LOPES trace (i.e. f re f = 0), which is deter- 
mined by the KASCADE-Grande trigger, because it starts the 
LOPES read out. As it simultaneously triggers the pulse gen- 
erator of the delay measurement, the condition f re f - to = const 
is fulfilled, and the delay f can by obtained as the arrival time 
t t of the calibration pulse in the trace of the calibration event: 
f = t t - = t t . 

This pulse arrival time f t is determined in a subsequent ana- 
lysis as time of the positive maximum of the up-sampled trace 
(like shown in fig. |5J- When repeating the measurement for 
the same channel several times, the measured pulse arrival time 
f t is stable within about one sample of the up-sampled trace 
(RMS of 12 successive events < 0.1ns, trace up-sampled to 
(12.5/2 7 w 0.1) ns sample spacing), if the amplitude of the cal- 
ibration pulse is chosen high enough for a sufficient signal-to- 
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Figure 5: Example for the delay measurement: the relative delay between two 
antennas (solid blue and dashed red) can be measured as the difference of the 
times when a calibration pulse is received, which has a height of about 1 V and 
is fed directly into the antenna cables. The relative delay is mainly caused by 
different cable lengths. 



noise ratio. Hence, this measurement method enables us to de- 
termine the relative delays Ar m n with a statistical error of about 
V2 ■ 0.1 ns* 0.15 ns. 

Furthermore, systematic errors of the delay measurements 
have been studied in several ways, e.g., by repeating the mea- 
surements. Measurements of the relative delays performed on 
two consecutive days deviate by (0.4 ± 0.3) ns from each other 
(mean and standard deviation of 10 measurements). As another 
check for systematic effects, the pulse arrival time f t have been 
determined in four different ways, namely as time of the posi- 
tive maximum of the trace, the negative maximum of the trace, 
the maximum of a Hilbert envelope of the trace, and the cross- 
ing of half height of a Hilbert envelope of the trace. The statisti- 
cal error of the relative delays is about the same for each method 
(~ 0.15 ns). But the value of the relative delays AT m41 depends 
on the way the pulse arrival times f t are calculated. Only the rel- 
ative delays calculated by the positive and negative maximum 
of the trace agree within the statistical error of about 0.15 ns. 
The relative delays calculated by the maximum of the envelope 
and the crossing of half height of the envelope disagree slightly 
with each other, and the delays calculated by the positive or 
negative maximum of the trace are highly inconsistent with the 
delays calculated by the maximum of the envelope, as they all 
have a statistical error of about 0.15 ns, but differ by up to a few 
nanoseconds. 

Under the assumption that the electronics of all channels be- 
haves identically, all methods for the determination of the pulse 
arrival times should lead to exactly the same relative delays. 
Hence, the explanation for the observed inconsistency is that 
the properties of the different channels are not exactly the same. 
Indeed, after correction for all measured differences, namely the 
amplification factor and the dispersion (see next section), the 
inconsistency between the delays obtained from the different 
methods is reduced. But still, there remains a deviation of up 
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Figure 6: Deviations between delays f calculated by the negative maximum 
of the up-sampled trace and the maximum of the envelope, mean shifted to 0, 
standard deviation =1.7 ns. The histogram contains 30 deviations of one delay 
measurement campaign of all 30 LOPES antennas. 



to a few nanoseconds for some channels, and the average devi- 
ation between the relative delays calculated by the maxima of 
the trace and by the envelope of the trace is of about 1 .7 ns (see 
fig. |6j. This shows the difficulty to fully correct for different 
channel properties. Or in other words, in designing the elec- 
tronics for a new radio antenna array one has to pay attention 
that components are from the same batches, etc. 

In the standard analysis of the shower reconstruction a cross- 
correlation beam is formed using the trace and not its envelope. 
Therefore, we have decided to use the delays calculated by the 
time of the positive or negative maximum of the trace of the 
calibration pulse. Thus, we minimize the systematic uncertain- 
ties introduced by the effect mentioned above. But still, there 
is another source of systematic uncertainty: The distance of the 
positive and negative maximum of the trace is about 9 ns be- 
cause the response of the bandpass filter causes an oscillation 
with the center frequency of the used band. This oscillation 
and, thus, the distance between positive and negative maximum 
and the resulting relative delays depend only little (~ 0.5 ns) on 
the shape of the calibration pulse. This could translate into 
a systematic uncertainty in the same order when determining 
pulse arrival times, if the pulse shape of cosmic ray radio pulses 
changes with lateral distance, as it is predicted by simulations 
I2l- 

Another check for systematic errors was to shift the emis- 
sion time of the calibration pulse by integer and non-integer 
multiples of the sampling clock, and no effect on the relative 
delays has been observed. This proves that up-sampling works 
reliably, and that the determination of the arrival time of radio 
pulses does not depend on how these pulses arrive relative to 
the original sampling clock. Consequently, neither up-sampling 
nor the original sampling rate of 12.5 ns do introduce any sig- 
nificant systematic errors. 

Summarizing, the total error on the relative delays AT m41 is 
below 1 ns for the standard cross-correlation beam analysis, 
which is more than sufficient for interferometric measurements 
with LOPES. For other analysis methods, like a lateral distribu- 



tion of pulse arrival times, the total error will be higher, due to 
the inconsistency of the different ways of calculating the pulse 
arrival time. In such a case the uncertainty is estimated to be in 
the order of 2 ns. 

The relative delays obtained by the described method are 
consistent with those determined earlier by solar burst mea- 
surements. The new method, however, has two fundamental 
advantages compared to using astronomical sources: The re- 
sulting delays do not contain any systematic uncertainty related 
to the errors of the measurement of the antenna positions, and 
the delay calibration can be done at any time. For LOPES, this 
is especially important because, due to the high noise level in 
Karlsruhe, solar bursts are the only astronomical source visible, 
and thus, continuously emitting astronomical sources are not 
available for calibration. The described method for delay mea- 
surement is repeated roughly once per year or whenever any 
changes in the experimental setup require it. 

5. Dispersion 

Dispersion is the frequency dependence of the group velocity, 
respectively of the group delay of a system. In case of dis- 
persion, waves at different frequencies propagate with different 
speeds, leading to a linear distortion of broad band radio pulses. 
For LOPES the dispersion of the analog electronics (which is 
mainly caused by the band-pass filters) has been measured with 
a network vector analyzer. Hence, the dispersion of the filter 
can be removed in the subsequent analysis by multiplying the 
appropriate phase corrections to the frequency spectrum of any 
recorded data. 

The effect of the dispersion has been studied with test pulses 
from a pulse generator which has been connected to the analog 
electronics instead of the antenna, like for the delay measure- 
ments. Different shapes of test pulses have been examined, and 
one example is shown in figure Qbefore (left) and after correc- 
tion for the dispersion (right). For most pulse shapes the disper- 
sion leads to a change in amplitude and FWHM of an Hilbert 
envelope of the up-sampled field strength trace of about 10 %. 
As the influence of the filter dispersion is largest close to the 
edges of the frequency band, the mentioned distortion effects 
can be reduced from about ten to a few percent, when using 
the sub-band from 43 to 74 MHz, only. For radio experiments 
with unknown dispersion such a selection of an inner sub-band 
would be a possibility to reduce systematic uncertainties origi- 
nating from pulse distortions. 

Because the radio pulses from real cosmic ray events are sim- 
ilar to the used test pulses (at least within the used frequency 
band), distortion effects in the same order of magnitude are ex- 
pected for real events (i.e. changes of a few percent of amplitude 
and FWHM). In addition the pulse arrival time changes by up 
to a few nanoseconds, depending on how it is calculated (e.g., 
value at pulse maximum or at the crossing of half height). These 
are changes of the absolute value which have a similar effect for 
all channels, as equal electronics is used, and thus, the disper- 
sion of each channel is approximately the same, which has been 
verified by measurement. Under the assumption that the cos- 
mic ray radio pulse shape does not change much on the lateral 
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Figure 7: The pulse distortion by the dispersion of the band-pass filter is visible in the measurement of a short test pulse. Crosses indicate the sampled data points. 
The lines correspond to the up-sampled signal and a Hilbert envelope of the up-sampled signal. In the raw data (left), the pulse is partially delayed by more than 
100 ns. After correcting for the dispersion in the analysis (right), the pulse becomes symmetrical and its height and width change by about 10 %. 



extension of LOPES (~ 200 m), it should be distorted by every 
antenna and its corresponding electronics in the same way. This 
means that the impact of the dispersion on the relative timing is 
expected to be much smaller than the observed absolute shifts 
of a few nanoseconds. Consequently, the dispersion of LOPES, 
even if not totally corrected for, should not spoil the capability 
to achieve a relative timing accuracy of about 1 ns. 

As all LOPES antennas are from the same type, their disper- 
sion is expected to affect the relative timing between the indi- 
vidual antennas only marginally. By this, it is acceptable that 
the dispersion of the LOPES antenna type is not known. It is 
difficult to measure, because the LOPES antenna can be used 
as receiver, only, and thus the two antenna method which is 
normally used for the determination of the dispersion, cannot 
be applied. In figure [2] the traces of a real cosmic ray event 
are corrected for the dispersion of the filters, and the remaining 
pulse distortion seems to be smaller than those shown in fig- 
ure |2l where the calibration pulse is affected by the dispersion 
of the filters, only. Thus, the sum of the dispersion of all other 
components, including the antenna, is assumed to be lower than 
the dispersion of the filters. Nevertheless, due to the high noise 
level for real events, and because the exact shape of the cosmic 
ray radio pulses is unknown, this can not be expressed quanti- 
tatively. 

For LOPES STAR which uses different antennas and electron- 
ics, the dispersion of the complete system has been measured 
1 2 ill . It was found, that the dispersion of the cables can be 
neglected, but the dispersion of the antenna itself cannot. It 
can be of the same order of magnitude as the filter dispersion. 
For this reason, future radio experiments should aim either for 
antennas with low dispersion or for antennas with well-known 
and, thus, correctable dispersion. Correcting pulse distortions 
induced by the antenna dispersion is especially important for 
larger-scale antenna arrays, if it turns out that the cosmic ray 
radio pulse shape changes with lateral distance. This could also 
have implications for the application of interferometric analysis 
methods, e.g., forming a cross-correlation beam. Hence, larger 




Figure 8: The dipole antenna used as beacon to monitor the timing of LOPES. 



antenna arrays, like AERA at the Pierre Auger Observatory 
or LOFAR [8], have the opportunity to test this. 



6. Monitoring the timing with a beacon 

Experience with LOPES has shown that the timing is not ab- 
solutely stable. Instead, once in a while, jumps by one or two 
clock cycles (12.5 ns) occur. In addition, small drifts or changes 
of the relative delays, e.g., with changing environmental tem- 
perature, cannot be excluded, as the electronics has not been 
designed for sub-nanosecond stability. Independent of the rea- 
sons, any changes of the timing have to be accounted for, to 
achieve an overall timing accuracy in the order of 1 ns. As the 
exact variations of the delay are not predictable, a continuous 
monitoring of the timing is needed which provides the ability 
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to correct the timing in the subsequent analysis on an event-by- 
event basis. 

For this monitoring we have deployed an emitting dipole an- 
tenna, a beacon, on top of a building of the Karlsruhe Insti- 
tute of Technology, at about 400 m distance to the center of 
LOPES (figs.[8]and[T]). This beacon permanently transmits two 
sine waves at constant frequencies of 63.5 MHz and 68.1 MHz 
(width < 100 Hz) at a low power of -21 dBm (» 0.008 mW). 
Thus, every LOPES event contains a measurement of the phases 
at these frequencies, which can be obtained by a Fourier trans- 
form into the frequency domain. Any variation in the relative 
timing between two antennas can be detected as a variation of 
phase differences at each beacon frequency. 

The phase of the continuous beacon signal at an antenna 
depends on the distance and the orientation angle of the an- 
tenna towards the beacon as well as on the delay of the cor- 
responding channel. Let us assume for a moment, that there 
are two antennas at an equal distance and angle to the beacon 
and with an equal delay. If we consider just one beacon fre- 
quency, e.g. 63.5 MHz, the two antennas would measure the 
same phase at this frequency (except for small deviations due 
to noise). Thus, a variation in the relative delay between the 
two antennas would immediately lead to a change of the mea- 
sured phases. If, e.g., the relative delay shifts by 1 ns, the differ- 
ence between the measured phases at the two antennas would 
be A(f> = Ins -63.5MHz -360° = 22.9°. Correspondingly, a 
measured phase difference can be converted in a shift of the 
relative delay. 

Now let us consider the more realistic case, that we have two 
antennas with different angle and distance towards the beacon 
and different electronics and cable delays. As the distance and 
the effect of the antenna orientation is not precisely known (be- 
cause there is no need), we expect to measure a different phase 
at both antennas at the beacon frequency. As long as neither the 
distance, nor the orientation, nor the relative delay do change, 
the difference between the phases measured at both antennas 
A0, e f would be arbitrary, but constant. Thus again, changes of 
the relative delay can be detected as changes in the phase dif- 
ference A<p. The only difference to the case above is, that these 
changes of the phase differences will happen not with respect 
to 0°, but with respect to A0 re f. 

The important point is to define A0 re f for each antenna with 
respect to a fixed antenna as reference (arbitrary choice) and 
each beacon frequency when the delay is exactly known. There- 
fore, we determine A<p le f as an average of the events taken at the 
time when we do the delay calibration described in section [4] 
This way we can monitor and subsequently correct any varia- 
tion in the timing back to the values obtained in the delay cali- 
bration. 

The limitation of the accuracy of the measurement of the 
phase differences is given by the noise and by systematic ef- 
fects. The noise of the phase measurement depends (within 
reasonable limits) on the signal-to-noise ratio of the beacon sig- 
nal, where the amplitude of the beacon emission can be chosen 
such that a sufficient accuracy is achieved. In case of LOPES 
we have chosen to emit each frequency at -21 dBm. The noise 
of the phase measurement has been determined by the jitter of 
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Figure 9: Phase differences between two antennas at both beacon frequencies 
for the first ten events per day for a whole year (May 2008 - May 2009), ex- 
cluding a few days of down time. The absolute value of the phase differences 
is not meaningful. Though, the changes are, which amount at both frequencies 
to about 1.5 ns (~ 35°) between summer and winter. Details see text. 



the phase differences in successive events and corresponds to an 
accuracy in the order of ~ 0.3 ns. Aside from that, the additional 
noise introduced by the beacon signal to the data is negligible, 
as the cosmic ray radio pulses are broad band and extend over 
the entire frequency spectrum. On the other hand the beacon 
signal is visible only in a few fixed and defined frequency bins 
and can be suppressed by artificially reducing the amplitude at 
these bins in the data analysis or in the hardware of the trigger 
logic if a radio self-trigger system is applied. 

In figure |9j the phase differences at both beacon frequencies 
between two LOPES antennas are shown for the first ten events 
of each day for one year. An annual drift of the phase differ- 
ences which corresponds to about 1.5 ns (~ 35°) can be seen 
consistently at both frequencies. The reason for this annual 
drift is not definitely known, but might be due to environmen- 
tal effects, in particular changing temperature, as the effect is 
largest in summer and winter. Also a jump in the timing of 
two clock cycles (25 ns) is visible which occurs during one day. 
Here it becomes obvious, that at least two beacon frequencies 
are needed, as changes in the timing larger than half a period 
(~ 9 ns) could otherwise not be detected unambiguously. A con- 
sistency check between the results at both frequencies is also 
necessary to identify a few noisy events (like the outlier in the 
bottom left corner of figure|9]), for which the beacon correction 
of the timing cannot be performed. 

When inspecting the data carefully, some features in the plot 
of the phase differences can be seen which do not occur simul- 
taneously at both frequencies - contrary to the visible general 
drift. These features are due to systematic effects and in prin- 
ciple decrease the achievable timing accuracy. Possible reasons 
for systematic effects are changes in the emitted beacon sig- 
nal (e.g., if the frequency generation is not absolutely stable), 
changes in the propagation of the signal from the beacon to the 
LOPES antennas (e.g., due to different atmospheric or ground 
properties), and non-random (e.g., human-made) noise at the 
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Table 1 : Summary of uncertainties of the relative timing and their relevance for the interferometric cross-correlation beam (CC-beam) analysis 



Effect 


amount [ns] 


reducible by 


to [ns] fl 


relevant for CC-beam 


ADC sampling frequency 


12.5 


up-sampling 


<o.i 


no 


measurement of relative delays: 

• repetition on consecutive days 

• different shapes of calibration pulse 

• method of pulse time determination 


0A h 

0.5 

1.7 




0.4 

0.5 
1.7 


yes 
parti ally c 
no 


remaining dispersion after correction 
for known properties of the filter 


< 1" 




< 1 


partially'" 


variations of the timing with time: 

• occasional clock jumps 

• drifts (e.g., summer vs. winter) 


1 or 2 samples 
(12.5 or 25 ns) 
up to 1 .5 


beacon 
beacon 




~ 0.7'' 


(yes) 
yes 


total uncertainty (quadratic sum) 

total uncertainty for cross-correlation beam 


up to ~ 28 




~ 2.0 
0.8-1.4'' 


no 

yes 


"The uncertainties due to the delay measurements and due to dispersion are 
'contains statistical error of 0.15 ns. 


not further reduced. 







c depends on how much the shape of the air shower radio pulse changes with lateral distance. 

d The exact amount is unknown, but assumed to be significantly smaller than a few nanoseconds (see section|5J- 

""value determined from cross-check of beacon correction with delay measurements (see section[6): 

The main contribution is a deviation of 0.6 ns between the results of both beacon frequencies which contains a statistical error of 0.3 ns for each frequency. 



beacon frequencies. As the scale of the observed features is 
significantly smaller than 1 ns, they have not been investigated 
in detail, and should not limit the ability of the beacon method 
to achieve a timing accuracy of 1 ns. 

As a cross-check, the changes of the delays between two 
dates roughly one year apart, have been measured with the 
method described in section [4] and compared to the changes of 
the beacon phase differences between the same two dates. The 
relative delays measured with the method of section [4] changed 
by (0.6 + 0.4) ns between the two dates (mean and standard de- 
viation of the absolute change of all 30 antennas). This itself is 
not unexpected as the electronics was not designed to be stable 
on a sub-nanosecond level. Comparing these changes of the de- 
lays measured with the method of section |4] with the changes 
observed by the beacon, reveals some systematic effects. In 
the ideal case, the phase differences at both beacon frequen- 
cies should change by exactly the amount corresponding to the 
changes of the delays. In reality, the changes observed at the 
two beacon frequencies are not totally equal, but the changes 
of the phase differences at the first beacon frequency and the 
second beacon frequency differ by (0.6 ± 0.3) ns. This is larger 
than the statistical error which is about 0.3 ns (see above). The 
changes observed at both beacon frequencies have been aver- 
aged, to check, if they are consistent with the changes of the 
delays measured with the method of section |4] the changes de- 
termined by both methods deviate by (0.7 + 0.5) ns from each 
other (average of the individual deviations of all antennas). 

Hence, systematic effects on the beacon signal seem to play 
a role, and it cannot be excluded that the observed drifts of the 
phase differences are - at least partly - not due to drifts of the 
electronics or cable delays, but due to these systematic effects. 
Nevertheless, this does not undermine the ability of the beacon 
method to monitor and correct real changes in the timing, like 
the described clock jumps, and to provide for each event a tim- 



ing accuracy in the order of 1 ns which is required for digital 
radio interferometry. 

Finally, a beacon cannot only be used to monitor the timing 
of an antenna array, but is valuable to check the health of the ex- 
perimental setup in general. As it provides a defined reference 
signal visible in each event, most possible failures of the an- 
tennas or the electronics are detectable by monitoring the bea- 
con signal. For example, we have been able to exactly find the 
date when we accidentally switched the cables of the two po- 
larization channels of one antenna, by investigating the phase 
differences at the beacon frequencies between these channels. 

7. Conclusion 

The methods described for the time calibration of LOPES are 
especially useful for radio antenna arrays in a noisy environ- 
ment, where the calibration with astronomical sources is not 
possible. They allow the determination of the electronics and 
cable delays with a very high precision, which can in principle 
be below 0.5 ns. Systematic effects, however, limit the actual 
achieved accuracy of the delay measurement to below 1 ns for 
our standard, interferometric cross-correlation beam analysis 
and to about 2 ns for the direct measurement of pulse arrival 
times. In addition, the dispersion of the electronics has been 
measured and is taken into account in the analysis of cosmic 
ray air shower radio pulses, to avoid systematic uncertainties in 
the pulse height which can be up to 10 %. 

Furthermore, we continuously monitor any variation of the 
timing with narrow band reference signals from a beacon, thus 
achieving an overall timing accuracy in the order of 1 ns for 
the cross-correlation beam analysis (see table[T]i. This way the 
nanosecond time resolution required for digital radio interfer- 
ometry is achieved for each event, and the phased antenna array 
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LOPES can be used as a digital interferometer which is sensi- 
tive to the coherence of the air shower radio emission. 

Finally, monitoring of the timing with a beacon is an interest- 
ing feature for any radio antenna array. As in principle the phase 
differences at the beacon frequencies are sensitive to any varia- 
tion of the relative timing, even the timing accuracy of antenna 
arrays without stable clocks should be improvable to about 1 ns. 
Hence, a beacon should provide any radio experiment in the 
MHz regime with the capability to do interferometric measure- 
ments. For example, the application of a beacon and the pos- 
sibility of interferometric measurements of the cosmic ray air 
shower radio pulses with larger arrays is presently investigated 
at the newly developed antenna array AERA at the Pierre Auger 
Observatory. Also LOFAR will apply the described methods 
for time calibration, and observe the radio emission of cosmic 
ray air showers with a much denser array. 
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